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Figure 1: Sketch of the model. 


Protostellar jets and magnetic diffusion 

Miljenko Cemeljic^ and Christian Fendt^’^ 

^Astrophysikalisches Institut Potsdam (Germany), ^Universitat 
Potsdam, Institut fiir Physik; email: cemeljic,cfendt@aip.de 

Abstract 

We investigate the evolution of a disk wind into a collimated jet under 
the influence of magnetic diffusivity. Using the ZEUS-3D code in the 
axisymmetry option we solve the time-dependent resistive MHD equations 
for a model setup of a central star surrounded by an accretion disk. The 
disk is taken as a time-independent boundary condition for the mass flow 
rate and the magnetic flux distribution. We find that the diffusive jets 
propagate slower into the ambient medium, most probably due to the 
lower mass flow rate in axial direction. Close to the star we find that a 
quasi stationary state evolves after several hundreds (weak diffusion) or 
thousands (strong diffusion) of disk rotations. Magnetic diffusivity affects 
the protostellar jet structure de-collimating it. We explain these effects 
in the framework of the Lorentz force. 
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1 Magnetic jets from accretion disks 

Astrophysical jets detected so far seem to be attached to objects 
where an accretion disk is indicated to be present. The similarities 
between jets from the different sources imply that the basic mecha¬ 
nism for jet formation should be the same. 

What kind of mechanism turns the inflowing matter of the ac¬ 
cretion disk into an outflow from the disk (or the star) is still not 
really known, although it seems to be clear that magnetic fields play 
a major role. 
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Our model is sketched in Figure |^. In order to model the time- 
dependent evolution of jet formation, we solve the set of resistive 
MHD equations, 
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We apply a polytropic equation of state, p — Kp^ with a poly¬ 
tropic index 7 = 5/3. We do not solve the energy equation but 
define the internal energy of the system with eq. (4). We may ex¬ 
pect that the turbulence pattern in the disk may also enter the disk 
corona and the jet, and that the jet flow itself is subject to turbulent 
diffusion. The magnetic diffusivity is denoted by the variable p, and 
other have its usual meaning, with pa = p//3t with Pt = const., an 
Alfvenic turbulent pressure, added to the hydrostatic pressure p. 

2 Results and discussion 

The results of one reference set of global simulations with high res¬ 
olution (numerical mesh of 900 x 200 grid points, physical grid of 
{z X r) = (280 X 40)ri) are presented in Figure 
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Figure 2: Global evolution of the jet on a grid of (zxr) = (280x40)ri with a 
resolution of 900x200 elements. Shown is the state of evolution after t=400 ro¬ 
tations of the disk inner radius for different magnetic diffusivity, 77=0,0.01,0.1 for 
top, middle and bottom panel, respectively. Colors (blue to yellow in decreas¬ 
ing manner) indicate density. Lines denote twenty linearly spaced magnetic flux 
surfaces (or poloidal field lines). The figure demonstrates that the bow shock 
advances slower with increasing diffusivity. 
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Figure 3: The evolution of the inner jet approaches the quasi-stationary state. 
Shown are poloidal magnetic field lines in the case of rj = 0.1 for different 
time steps, t = 250, 300, 350, 400 (thick solid, thin solid, dashed and dotted 
lines, respectively). In these simulations beta plasma is half the value from 
the global simulations above, but this just changes the position of the Alfven 
surface above the disk slightly. Grid size is 280x80 elements for a physical size 
of (140x40)ri. The picture shows how the poloidal magnetic field lines diffuse 
outwards but approach a (quasi)-stationary state after 400 rotations (see the 
dashed and dotted lines almost coinciding). 

In order to investigate only the gross behavior of the jet flow and 
not its structnre in detail, we rnn another set of simulations with 
lower resolution (numerical mesh of 280x80 grid points, physical 
grid of (zxr)=(140x40)ri) - see Figure 


In apparent contrast with the bow shock propagation from Fig¬ 
ure H is the increase of the jet velocity with diffusivity (Figure H, 
left panel). 

Below the ‘’critical rj‘’ in our choice of parameters (77 < 0.5), 
the time of reaching the quasi-stationary state of the inner jet is 
approximately linearly dependent on rj (Figure^, right panel). 

Interesting is to note that de-collimation of the poloidal mag¬ 
netic field lines is not followed by the poloidal velocity vectors - see 
Figure ^ 

The much faster computation of the low resolution simulations 
allowed us to follow the jet evolution for a very long time even in the 
case of a high magnetic diffusivity (up to 4000 disk rotations). The 
observed de-collimation of the matter flow with increasing diffusivity 
is most evident if we plot the mass and momentum fluxes across the 
boundaries of the inner jet region. 

In combination with various physical effects - magnetic and iner¬ 
tial forces, pressure and gravity - the magnetic diffusion will modify 
the MHD structure of the jet. As the flow approaches the Alfven 
surface, a toroidal magnetic field component is induced (‘’wound- 
up‘’) due to the inertial back-reaction of the matter on the field. 
The toroidal field may lead to (de-) accelerating Lorentz forces 
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Figure 4: Left: Components of the velocities. Cuts in Z direction at r=15 for 
the quasistationary time of simulations. Poloidal components are given in solid 
lines, toroidal in dashed lines. Two different ordinate scales are used to expose 
relative changes of the velocities. Left ordinate gives the values of the poloidal 
velocity, when right gives the values of the toroidal velocity. Both components 
increase with the increasing diffusivity. Right panel: The times when the quasi¬ 
stationary state is reached, in dependence on diffusivity. 
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Figure 5: The de-collimation of the stationary state poloidal magnetic field due 
to magnetic diffusivity. Thick black lines for 77 = 0 , and thin lines for 77 = 0 . 1 . 
Overly of the poloidal velocity vectors show that the velocity pattern is not 
de-collimated. Thick blue line is the Alfven surface, and thin blue line is the 
fast magnetosonic surface for the diffusive simulation. 
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Figure 6: Time evolution of the mass flow ratio between the radial out¬ 
flow boundary and the axial outflow boundary in the inner part of the jet, 
(zxr)=(60x20)ri. Parameter ?7=0,0.1,0.5 in solid, dashed and dotted line, re¬ 
spectively. For higher diffusivity, the mass flux ratio in the quasi-stationary 
state increases, indicating a decrease in degree of collimation. 


-^L,|| ~ j± X and (de-) collimating forces ~ j\\ x B where, 
here, the perpendicular and parallel projection is made with respect 
to the poloidal magnetic held (which, only in the case of ideal sta¬ 
tionary MHD is parallel to the poloidal velocity). 

The winding-up of the poloidal magnetic helds is less efficient 
since magnetic diffusion leads to a slip of matter across the held. The 
induced toroidal magnetic held is weaker leading to a less efficient 
acceleration by Lorentz forces but also to a de-collimation. As a de- 
collimation of the poloidal magnetic structure also implies a smaller 
launching angle for the sub-Alfvenic how, the magneto-centrifugal 
acceleration mechanism may work more effective. 
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Figure 7: Lorentz forces in the jet for different magnetic diffusivity 77 = 0, 0.1, 0.5 
(solid, dashed and dotted lines, respectively). (Normalized) values of the force 
component perpendicular (left panel) and parallel (right panel) to the field line. 
The toroidal component is similar to the parallel one. Components are projected 
along a flux surface leaving the box of the inner jet close to (R=20,Z=60)-corner. 
The sign is defined as follows. For the perpendicular component the positive 
sign denotes the r-direction (de-collimating force). For the parallel component 
the positive sign denotes the z-direction (accelerating force). 
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